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Prof. A. TANAKADATE (1857-1952) 


On the morning of May 16, 1952, Prof. A. Tanakadate, the honorary member 
of our society, was going to attend an semiannual meeting of the Acoustical Society 
of Japan, but feeling his left hand lightly numb, he took a rest at home. A symptom 
of brain disease (cerebromalicia) appeared gradually, accompanied by unconscious 
sleep for three days, and towards the evening on 21st he quietly closed his century- 
long life, which was a lofty and faithful sequence of devotions to human culture. 

Prof. Tanakadate was born on Sept. 18, 1857’ at Hukuoka in Iwate prefecture, 
a small castle-town near the north end of Honsyu-Island. It was just the beginning 
of the era of emancipation of Japanese people from a long’ continued policy of 
seclusion from the world. On entering Tokyo University in 1878, he gave up his 
young ambitious hope to become a statesman, and decided to learn physics, in which 
he recognized the essential nature of the forthcoming civilization. In the course 
_ of his study his inborn qualities for a scientist were developed by excellent teachers, 
_ Prof. Thomas C. Mendenhall and Sir Thomas Alfred: Ewing in Tokyo University, and 
_ Sir William Thomson (Lord Kelvin) in the University of Glasgow (1888-1890). From 

1891 to 1917, he kept a chair of professorship in physics at Tokyo University. A number 
of able students who learned physical science under his severe training laid foun- 
dations to scientific researches of the present Japan. He was one of the founders of 
the Earthquake Investigation Committee (1892), the Mizusawa Latitude Observatory 
(1899) and various geophysical organizations. As early as 1910 he succeeded in taking 
_ instantaneous photographs of the whirlwinds generated by airscrews in his first wind- 


tunnel, and afterwards stimulated the establishment of the Aeronautical Research Insti- 
tute, which was realised in 1918 and originated the fundamental researches for the 


aeronautics and the airplane technology in Japan. He retired from the professorship - 
in 1917, but his life as a leader of science became more active and important than 
r. He devoted his efforts to the spread and circulation of scientific knowledge over 
1e world as well as within the country. A considerable amount of his energy was 
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which was often very useful to soften the occasional tensions between peoples in the 
international congressess aS well as in the meetings at home. 

His enthusiasm for the international collaboration was not exhausted even in 
the later days. The writer received a letter from him just before starting from 
Japan for Brussels last year (1951) to attend the 9th General Meeting of IUGG. Some 
passages are quoted here: 

“The very city Brussels recalls to my mind the matorset cael event, the 

Inaugural Meeting of the International Union of Geophysics and Geodesy in 

1919. Most of the magneticians of my intimate acquaintance, such as Schuster, 

Chree, Bauer, etc., with whom I had the pleasure to collaborate offering my 

meagre efforts in starting the International Association of Terrestrial Magnet- 

ism and Atmospheric Electricity (then so called), are now regrettably passed 
away.” 
(then he describes two events in this meeting; the discussion about the name of 
IUGG, and the vote for the Exclusion Pact against Central Powers, closing with the 
following.) 

“Happy to say, such cloudy concepts are now clearey gone away with, and the 

true aspect of International Union is seen gleaming its future brilliance in 

promoting the progress of human knowledge. I wish to congratulate in this 
sense the Brussels Meeting 1951 as a surviving member of the Union. Please 
don’t forget to remember me kindly to my esteemed colleagues, Fleming, 

Chapmen, Maurain etc etc when you have occasion to talk with them,...... i 
He was appointed first presiddent to the established association (1919), and remained 
to fulfil his duty in the Executive Committee of the association until 1950. 

Although, in the enlightening days of Japan, his activity, as the leader of leaders, 
covered a wide field of science and technology, his personal interest, it seemed to the 
writer, was specially directed to the researches on the terrestrial magnetism 
and electricity. In 1884 he devised a set of instruments of measuring D with electro- 
magnetic method (Proc. Phil. Soc. Glasgow, 1889). The Tanakadate magmetometer, 

a nice and ingenious instrument for magnetic survey, brought a marked step forwards 

in the accuracy and the time saving of the magnetic work in Japan. The first 
complete magnetic survey of Japan and the initial routine observation of magnetic 
field in Tokyo (1897) were as well carried out by using this magnetometer. His report 
entitled “‘ A magnetic survey of Japan reduced to 1895.0 and sea level” (Journ. Coll. 
Sci. Tokyo Univ. 1904) is highly esteemed as a text book of magnetic survey up to 
the present in Japan. In the last chapter of the report, prof. Tanakadate assumed 

a magnetic image which caused the observed regional anomalies of the magnetic 
field, and sought its position and intensity from observed data in Japan, England and 
Hungary. It is truly interesting that by means of such a direct method of deter-- 
mination he obtained something about. 2000 km as the value of depth of the under- — 
ground images in the case of “single magnetic poles.” 

He took an active charge in the project of the II Polar Year as an officer of 


5; 
4 
ee abet. tie 


‘ 
- 


Prof. A. TANAKADATE 43 


IATME, and as the chairman of the Japanese Committee of the II Polar Year. In 
his report on the sudden commencements of magnetic storms (Trans. Lisbonne 
Meeting, IATME), we find valuable hints, to be taken into account in the case of 
similar plans of international cooperation in future. Many systematic magnetic 
observations were undertaken for survey, at the time of solar eclipse, or for the 
purpose of earthquake investigations under his presidency during the period 1880- 
1930 in Japan. 

Imperial Academy nominated him a life member in 1906, and since 1917 he 
assumed the title of professor emeritus of Tokyo Imperial University. From 1925 to 
1946, he was elected three times successively member of the House of Peers as a 
representative of the Imperial Academy and contributed so much to the cultural 
policies of Japan. Many honours given him inland and abroad celebrated his scientific 
career which culminated in his acceptance of Bunka- Kunsho, the highest honour 
granted to learned people in Japan. . 

Contrary to his severe and rigorous criticism for scientific works, he was 
endowed with a character of rare warmheartedness ever seen. His married life 
was short, lasting only one year in 1894 until his wife died, leaving him a daughter 
alone. Since then he did not marry again. But he was never lonely; he was helped 
by a self-sacrificing partner, his daughter, Mrs. Mineko Tanakadate at home, and a 
great deal of his love was poured to science, scientific friends home and abroad and 
his research successors. Until before the last illness he had never neglected to be 
present at meetings of various scientific societies. Everywhere we found him sur- 
rounded by a welcoming eiicd of scientists, young and old. 

By the death of Prof. A. Tanakadate we lost the father of ee science of 
geomagnetism and geoelectricity in Japan and a man of the noblest spirit, whose 
long life was wholly devoted to the love of humanity and search for the truth of 
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The Electron Density of the D Region 


By Teruo SATO 


Geophysical Institute, Kyoto University 


Abstract 


In this paper we study the electron density, its seasonal and daily 
variation, and the distribution for the height of the D region. For 
this purpose, the seasonal and daily variation of the reflection coe- 
fficient of the waves of 16kc/s and 150/s incident in the ionosphere 
were calculated assuming the distribution of the electron density, and 
maximum electron density N, (though its existence is unceatain) of 
the D region and the collisional frequency » at the height of N>. 

To calculate these values, the levels of the reflection of these 
waves were firstly discussed. Then No, », and the distribution of the 
electron density which give the reflection coefficient consistent with the 
results observed not only by the long wave but also by the short 
wave were deduced. It was found that 
(1) waves were reflected from the level x=sin*6[y,?+Z7]” and not 

x=sin?0 or x=sin?0-Z, 
(2) the distribution of the electron density has the form such as 
gy Mbcianen'arc ”) or a similar one, 

' (3) the electron density at the height of 90km and 75km is 1.8% 104 

/cec. and 6x10°/c.c. respectively in daytime in the summer and 
9x 10%/c.c. and 3x10°/c.c. in daytime in the winter, 
(4) if the region has the peak it is situated at the height of 92km- 
95km in daytime, - 
(5) and the collisional faequency at 92km is 4x 10°/s. 
§ 1. Introduction 

It has been observed that there are two types of the D regioh. One is the 
regular region situated under the E layer at the height of 70km-95km. The other . 
is a sporadic region which is situated at the height of 50km-60km, and which often 
reflects partially the wave incident in the upper layer. Ellyett“ (1947) summerises 
reports for the latter region observed by many workers. In this paper, we study 
about the former region. 

Many observations were made by several workers. [J. Holligworth® (1926), R 
Naismith (1931), S. Namba“ (1931), J.E. Best et al. (1936), K.G. Budden ef al. 
(1939), J.N. Brown e¢ al.” (1950), A-H. Benner® (1951), and R.A. Helliwell e¢ a1. 
(1951)]. Frequencies used by them are in a range of 16kc/s-325kc/s, and the height 
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where these waves are reflected is above 7Okm. In these works it had not been 
determined whether the D region is the tail of the E layer or not, and whether the 
maximum electron density of the D region exist or not. It was the main aim for them 
to observe the height of reflection, the polarization, and the reflection coefficient of 
the wave. 

The reflection height and the reflection coefficient or long waves at the frequencies 


=~ / oF ice /a « ; ; 
16kc/s-325kc/s are summerised in Table 1, and the deduced electron density at the 
height of 70km and 90km are shown in Table 2. 


Table 1 
| angle of | hight of refletion | 
wees frequency aaniest (km) reflection 
(kc/s) | ; period 
(degree) day night coefficient 
Benner 150 0 ¥ | 0.05 (midday) 
“7 99~98 | 0.5 (midnight)| 1939-2—1950.2 
: | "325 e 
Helliwell } 0 90~130 
100 a 
et al. he 0 Re 0.88 1947.2—1948.8 
Brown Watt | 50 0 80~76 80 1950.1 
ene 72~75 
c ong distance (summer 
Hollingworth| 25 ae 1924. 6-1925-12 
(autumn) 
Best, Budden, 1934.9 
ilies, 16 ere 75~95 92 | 0.25 (day) | 193).11 
Ratcliffe (Sep. ) (Nov.) 0.5 (might) 1936.8 
1937.9 
17.4 Long distance | 80 or 100 a 
Namba 19.8 Transmission (assume) 80 or 100 ; 1930.4 
Tabel 2 
Shes. S Worker 
Budden Bates (0 ; 3 
Height ~~ et al. Seaten Bremmer (1) | Mitra(@? 
70km Boe Luce aise: OPK LU 225 x 108 102 
90km 1.5103 2x 104 


It is difficult to make clear a question as for the existence of the maximum 
electron density and its magnitude in the present step of knowledges. Aim of our 
study is to deduce the distribution of the electron density and its value at any height 
of the D region which gives the observed reflection coefficieut and reflection height 
of long waves. In this study, it is convenient for the calculation to assume the max- 
imum electron density (N)) and the collisional frequency v, at the height of No. 

" If we get the value.of reflection coefficient consistent with the observed one 
using assumed N,, vp, distribution of collisional frequency and electron density, and also 
the observed results by the short wave are satisfied by these quantities, then these quan- 
tities are the required values and distribution. We used the formula in the calculation 
of the reflection coefficient which is derived by M.V. Wilkes (1940) and extended by 


W. Pfister‘ (1949). 
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§ 2. Calculation of Reflection Coefficieut | 
The equation for a plane wave incident in the ionosphere at an angle @ with 


the gronnd for a long wave is given by M.V. Wilkes (1940) as follows: 


AL(s) 
— ds? 


+k'[u*(s)—cos*6]L(s) =0, ey 


where L is the amplitude of the wave, s is the coordinate of the height, is the refrac- 
tive index and nat, W. Pfister extends this as follows. 
The height of reflection is given by 
| 2°—1| =sin’0. (2) 
If the reflection occurs in a parabolic form with the center on the point of reflection, 


w’—1 is given by 


pwo—1=— ik? Ks? (3) 
Pep seas 
VK Es e (et 4 (4) 


where H is the scale height, z is the height in the Be of the scale height, w = 2nf is | 
a angular frequency of the mate and c is the velocity of the light. (1) is transformed 


as follows: 
aL (€) 1 (be 
Bilin 2B LEN aa es 
; ds? a(t 2 jh (S)=0, 
3 vhs __sin* 1 Laaliee 
with tao + |pi¢ 
Cas 2-AVIK nt me ae nt ae (530 


wave and the incident wave. If the field strength and the gradient with the height 
are steady, the formula for the reflection coefficient is given by 


phate 1 dee 
wi 30 ‘Rksind ds L(€)o_ 


If we use the parabolic cylindrical function, ¢ is given by 


Stes ne da: ik gore = *| 
~Rsind ds ~V n+2 ; 
me 4» Pee 

mag 2) ; a 


et ey 
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where Serie a = N is the electron density of the D region, y is the collisional 
frequency and e, m are the charge and mass of electron. In this case, it is easily 
known that both ordinary and extraordinary waves are reflected at the same level, 
and it is difficult to separate these waves definitely. 

M.V. Wilkes“ (1947) assumes that a very long wave is not reflected at the level 
x=1, but reflected at the level x=1+yz, where yn 2H, Hy; is a component of the 
earth’s magnetic field vector along the propagation of the wave, and estimates that 
the electron density at the level of 75km is about 300/c.c.. RJ. Nertney“” (1951) shows 
that the wave is reflected at the level x=1 and x=1+yz (neglecting the collisional 
frequency), as the explanation of the stratification of the lower ionosphere which is 
observed by R.A. Helliwell e¢ al. (1951). Thus it is difficult to determine which of 
three cases is the true level of the reflection. We, at first, must determine a level of 
reflection, since it is necessary in the calculation of the reflection coefficient. At the 
level under the £ layer, the collisional frequency is considerably great, so that we 
cannot neglect the collisional frequency. In this case x=1-++-y,; must be corrected as 
x=[(+yz)?+Z]/*. For a long wave, yz>1, then x=[y,2+Z?]’. Thus the effect of the 
difference between the ordinary and extraordinary waves disappears and two waves 
are reflected at the same level. If we use an oblique incident wave, sin?@ must be 
multiplied to each term. 

In order to compare a calculated value with a measured value, we take two 
waves at the frequency of 16kc/s and 150kc/s respectively. Angle of incidence is 30° 


for the former and 90° for the latter. 


§ 4. Reflection at the Level x=sin*OLy2+Z°]/? 

Although the distribution of the electron density of the D region is not clear, 
it appears that this is a semilar distribution as that of upper ionosphere. Types which 
we presume are the Chapman’s, modified Chapman’s, parabolic, exponential and linear 
one. As we calculate the daily and seasonal variation of the reflection coefficient, it 
is desirable that the calculated value in any distribution of electron density does not 
much differ from that of the measured one. The daily variation of the reflection 
coefficient, when the parabolic, linear and exponential distribution are used, is abruot. 
So we use the chapman’s or modified Chapman’s distribution. 

If we assume, at first, as the electron density distribution of the D region, the 
Chapman’s type (and same for § 5 and § 6) which is given in the form 


1(1—2—e-*sexy) 
> 


N=WNoee 


and y=Yoe-* 


He Sree 
as the collisional frequency. Then L ++, | is given by 


ES Ep Sara ee eS las ee aoe ee 8 
|n+ z|=e sin [(e-#secy — 1)?+ 4e-*secz-sin’a]/? * b2+07)”" C8) 


ae 2 eed, 
(a) Case of 16kc/s err 
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If the earth’s magnetic field is 0.45 gauss, x=sin°6Lyz’ +Z]}”? is transformed as 


follows: 


= 4(1—z—e-*sexx) 


No. _—«O[5X 108+ vo%e=27}1/?* 


If we assume N, and v properly, the height of the reflection is deduced from (9). 


Substituting this value in (8), n+ is calculated. Thus, using (6), (7), we can cal- 


2 
culate the reflection coefficient. Lord for the wave of 16kc/s is given by 


= , sSina= 3 E Fits. 
[(e-*secy — 1)? + 4e-*secz-sin’a]/” [(voe™*)’+5 x BUN et Bi 
It is clear that the term involving v, in the denominator of the right-handed side of 


(9) is negligible of »)<10°/s, and is predominent if y)>10'/s. Therefore, it is not | 


1 | 3.464 ts yer? 


labourious work to deduce z when v< — 


i 

a 10°/s.. This is same for the calculation — 

ee of sina. 438%: { 
: eal The reflection coefficient for No= 
-< = 10'/c.c., 2x 104/c.c., 4x 104/c.c. and v= 
= : 2x10°/s 10°/s respectively, is given in 


Fig. 1. From the reason above men- 


73 tioned it is natural that the reflection — 
SEGX 
Fig. 1 Variation of reflection coefficient with 
secy for the wave of 16kc/c reflected J 
from the level x=sin?0{y,2+Z2}/2. siderably for the variation of secy. _ 
Angle of incidence is 30°. No is the (b) Case of 150kc/s. 5 
assume maximum electron density per 


‘cm, and vo is the collisional frequency 
. -per second at same height. 


coefficient does not vary much for the 
great variations of y) but varies con- 


jn+s| is given by — 


a5 | ie 3.76 x 20 © 
[(e-*secx —1)?+-4e-*secysin?a}/?* 


Ss 


N=lO* yal? 
Net erm 


LSNG2F U=5xI 
oi No= 2x10" ’ Yu=!0 
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The height of the reflection of the wave is derived from the equation 
eorliteee oe 
No ~ [Bx 10+ (e-#p,2}72 ° Se 


The reflection coefficient for N)=10*/c.c., 2x 10*/c.c.. 4X 10!/c.c. and »=10°/s, 5x 
10°/s, 10°/s is given in Fig. 2. Variation of the reflection coefficient for that of v» is 
small. In Fig. 2 the reflection coefficient of the wave for 150kc/s measured by A.H. 
Benner in January, which we compare with the calculated value, is shown. When we 
compare the curve in Fig. 2 and 3, with the measured curve it is easily found that 
the curve which we want cannot be gotten for above Np. In order to get a curve 
similar to the measured one, it is necessary to take a greater value of Ny and ». If 
we take N, as about 10°/c.c., it may be possible to get a necessary curve. But from 
a point of view that this value is almost equal to the maximum electron density 


in the normal E layer, we must try to get the coincidence in another way. 


§ 5. Reflection at the Level x=sin?0 
Collisional frequency is not involved in these calculations, so that, even if we 
determine the appropriate value of No and 1%, we would not be able to determine the 


height of the maximum electron density. 
(a) Case of oe 


For this wave ney is given by 
LT} 3.464 
ete erm pe 
The height of the reflection of the wave is found from the equation 
aie TEN 284, (14) 


By a similar procedure in the case of § 4, we find the reflection coefficient of the wave 
for N,=10?/c.c., 103/c.c. and 2x10%/c.c. in Fig. 4. It is found that the reflection 
coetficient in this case is remarkably 
great for small value of N,) compared 
with that in § 4. Since the electron 
density at the level x=sin76Ly Pua Ada “1s 
greater than that at the level x=1, the 
wave is to be absorbed more in the 
ionosphere in the former case than in 
the latter case. From the fact that the 
SEC x reflection coefficient in the former case 

Fig. 4 Variation of reflection coefficient with for greater value of ™ is ‘Smaller than 
ee ee that in the latter case for the smaller 


from the level «=sin®0. Angle of 
incidence is 30°. value of No, we suppose that the energy 


of the wave in the latter case is not Rates absorbed but apparently absorbed in the 
ionosphere, and thus the reflected energy is much less than’ that in the former case, 


S 


REFLECTION COEFFICIENT 
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Therefore, we cannot regard the level x=sin?@ as the level of the reflection of the wave. 


(b) Case of 150kc/s. 


For the wave | n to is given by 
4x9.41 
In ty 2 a “e-*secx—1 ° (15) 


The height of the reflection is found 


REFLECTION COEFFICIENT. 


from the relation 


1(1—2-—e-*secy) 


Noe ==. cal Op LG) 

; Fig. 5 Variation of reflection coefficient with 
Value for N?=10%/c.c., 5x10°/c.c., 10° secy for the wave of 150kc/s reflected 
/c.c. are shown in Fig. 5. from the level x=sin?6. Angle of 


For this wave, it appears that the incase ee 


suggestion of the apparent absorption of the energy of the wave Sie too. 
§ 6. Reflection at the Level x=sin0-Z 
By H.G. Booker, this holds only for the wave at the frequency of i i 


this case Jn+| is given by 


— > ae 


1 3.464 
eG Ok Be ak 12 OE 1 
n+ 2 | e-*secx+1 ° Alf 


The height of the reflection of the wave is found from the relation 


3(1—z—e-*sec x) 


Noe = 2.34 10-51. (18) aig y 


In Fig. 6 the reflection coefficient of the wave for N)=10#/c.c., 2x104/c.c., 4x10!/c.c. i 
and y¥)=10°/s, 10°/s, 104/s is shown respectively. We immediately find that the gradient Sa 


_ of the curve is contrary to the curves in Fig. 1-5. This implies that the absorption 


lation ek in § 4,§5 and” 
cee Bhat, it is appropraa, te 
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is necessary to take a large value of Ny and v, which we cannot allow. Thus we are 
forced to take another electron density distribution. According to the reason above 
mentioned, we use as an electron density distribution of the D region a modified Chap- 
mans type. As the gradient of the distribution must be slower than Chapman’s type, 
we assume that the distribution has for example the form such as N= Npe2 (1-2-2 *sec x) 
or a similar one. When we calculate the reflection coefficient for an assumed distri- 
bution for common Np-and yy) to two waves and get the results consistent with the 
results of the measurements for each wave, this Ny and vy» are the one which we 
require only in the standpoint of the long wave observation. 

IJ. Shaw“ (1951) show that the ionization density increases at a rate 1.71/km 
from the study of the cross-modulation of the wave. This is the value at the height 
of 85km-90km. On the other hand, D region is situated at the height of 75km-95km 
from the observations by many workers, and several workers suggest that the D 
region is not the tail of the E layer. So that the gradient of the ionization density 
given by Shaw may be that of the D region. When we select the form of an electron 
density distribution, we must note this fact. 


§ 8. Estimation of No and vo 


Now we come to the estimation of NM, and », using a modified form. Fig. 7 
nd Fig. 8 Show the results of calculation of reflection coefficient for two waves for 


S10 SSE aac = 
= ee 
= 08 S 
oe i 
wy iu 
~ 06 S 
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SEC 
Fig. 7 Variation of reflection coefficient with Fig. 8 Variation of reflection coefficient with 
secy for the wave of 16kc/s reflected secy for the wave of 150kc/s reflected 
from the level x=sin?6{y,?2+Z2}1/2. from the level x=sin20{y,2+Z2}1/2, 
Angle of incidence is 30°. Distribution Angle Of incidence is ZerO. ae 
of electron density is Nyes6-*-* “8%, of electron density is No#(-*— ™, 


differedt Nj and ». According to the measurement by J.E. Best ef al. for 16kc/s, the 
reflection coefficient is about 0.5 at night and 0.25 in daytime. Results of measure- 
ment by A.H. Benner is shown in F ig. 3. It is found that the modified distribution 
of electron density gives the value which we required. But we cannot yet determine 
N, and vo, since the same curves may be given for differedt MN, and vy», and distribution 
of electron density. 

It is convenient for the determination of NM) and vy to use two conditions. One 
is that the height of NM is at or above the height of 85km in daytime, and the other 
is the.supposition that the results of the measurement by short wave must be explained 
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by these quantities. F.R.W. White and L.B. Brown (1935) deduce the reflection coef- 
ficient of the wave reflected from the E, Fl and F2 layer and estimate the absorption 
by the D region. 

When we consider these conditions, our study is to determine Np and v ‘and 
distribution of electron density which satisfies following three conditions. 


Observed reflection coefficient (short wave)=f (No, vo) (a) 
Observed reflection coefficient (long wave) =g (No, vo, N distribution)  (b) 
Electron density, Reflection Height =h (No, vo, N distribution) (c) 


We start from (a), because this is the simplest way to determine unknowns. 

At first, we estimate the magnitude of N, and » by calculation of reflection 
coefficient using the formula of the absortion of the short wave. In order to avoid 
the complication of the integration, we use the Chapman’s type as the electron dis- 
tribution, for the D, Fl layer as well asthe E layer. Reflection coefficient for the 
ordinary wave is as follows; 


Are” Novocos? x ag) 


? 


Ron? me (oto) 


where suffixes d, e and f show the value concerning the D region, EF st F 1 beh 
respectively. 

When we derive the form (19), it is estimated that v=e-*. This form holds 
at the neighboughhood of any fixed level, for example, at the level of the maximum 
electron density of the D region and E and F1 layers. | 

The procedure is following. As the maximum electron density and its virtual 
height of the E and F1 layer are known, R, and Ry; are calculated assuming vo, a 
Yor and R, is calculated assuming Nog and voz for the wave reflected from the F2- 
layer. Same procedure is made for the waves reflected from the Fl and E aes ~ 


Table 3 


Te hv 
| Bei, ht Shaw“) Bates“ (20)| Best@) [Farmar®»)} Whites) way |Gibbon®| S.U.A. 
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Secondly, Ra+R.+Ryi or Rat+R, is compared with the value of the measurement 
respectively. In this process, it is: labourious work to select the correct collisional 
frequency by various height. Collisional frequencies at various height are shown in 
Table 3 and the example of this procedure is shown below. 


Example 


For the wave reflected from the F2 layer (4.5 mc/s) 


Total absorption (measured) 2.6 nepers Vo No 
Absorption by F1 layer 0.1 (2 x103/s) 2x 108 /e.c: 
Abshrption by F layer 12 (2.3 x 104/s) 1.45 x 10°/c.c. 
Absorption by D region 0.8 (4: x105/s) (104/c.c.) 
Sum Zall 

For the wave reflected from F1 layer (3.5mc/s) 
Total absorption (measured) 3.8 
Absorption by E layer 18 (2.3 x 104/s) 1.45 x10°/c.c. 
Absorption by D region 13 (4 x10°/s) (10°/c.c.) 
Sum a 

For the wave reflected from E layer (3.1mc/s) 
Total absorption (measured) 2:5 
Absorption by D region 1.5 (4 x10°/s) (104/c.c.) 


In this example, the collisional frequency given by K. Maeda is used and sum 
of the absorption by each layer is nearly equal to the value of the measurement 2.6, 
3.8 and 2.5 nepers respectively. From this example, the maximum electron density of 
the D region is about 10*/c.c.. But, we may take an arbitrary value of 4, so long as 
Rzt+ R.+Rei or Rz+R, is not greater than the measured value. Therefore we will 
get a various value of N, for the different yz. If we consider that the height of the 
maximum electron density is at or above 90km in daytime, »~¢ must not be greater 
than 4x10°/s when we use this distribution of the collisional frequency. Thus Nog 
will become greater than 10*/c.c.. For these Noa and voz, the reflection coefficient in 
Fig. 7 or 8 will become greater than the measured value so we must take another 
distribution of the electron density. By the repetition of this procedure, Nog and vog 
and eistribution of the electron density which satisfies (a) and (b) are gotten. 

But, on the other hand, the height of Nog is determined from the distribution 
of the electron density and the electron density at a certain height which is measured. 
If we take as v¢<4x10°/s, we are forced to take the sharper distribution than 
Noget (1 -2-e*8€CX) When this distribution is used, the height of NM) is- lower than 
the height corresponding the assumed ». If we select above quantities properly, (a) 
(b) and (c) are satisfied. 

Collisional frequency 4x10°/s is the value at the height of 90km, and this height 
is equal to the height deduced from both distribution Noge!-2-¢ "see and the 
electron density 2800/c.c. at the height of 85km which is derived from x=(y?+Z 7? 
and measument. Therefore, if we we use these values above three equations are 


54 T. Sato 


satisfied. This implies that this value is suitable to use and further the height of N, 
is 92km-95km. 

When the collisional frequency deduced from the measurement by other workers, 
the difference between the calculated and measured values become further great or 
small and there is the case that R,+Ry: or R, is greater than the measured value. 

F.G.W. Whitw and L.W. Brown shows that the absorption of the D region is 
1.15 neper in the winter. We can derive that Njz=9400/c.c., using cos?/?7 =0.193 and 
Yor=4X10/s. This is almost equal to the value calculated from our method, since the 
observation used here were made in daytime in September November and January. 

After all, we may show that Nog and voz are about 104/c.c. and 4-5 x 10°/s res ec- 
tively in the winter and at the same time, it is suitable touse the collisional frequency 
deduced by K. Maeda or D.R. Bates and F.W.S. Massy and to use the distribution of 
electron density expressed by Noe! (1—2—© *secx), 

§ 9. Seasonal Variation of the Electron Density 
It is presumed that the electron density of the D region increases in the summer 
similar to the E layer. By the measurement of A.H. Benner, the absorption in the 
winter is one-half of the absorption in the summer in unit of neper. As the increase 
of the absorption of the wave depends on the increase of the electron density and 
the collisional frequency, it is difficult to conclude that the increase of the absorption 
depends only on the increase of the electron density. According to the measurement 
by J.E. Best ef al., the level of the reflection is lower in the winter than in the 
summer. As magnitude of 3km is the order in a range of the error of the measure- 
ment and it is considered that the wave at the frequency of 16kc/s is reflected from 
the lower part of the D region, we may consider that the height of the maximum 
electron density does not change through a year. Then the magnitude of N) become 
about 2 x 10*/c.c. in daytime in the summer. This is deduced from the same procedure 
as that in § 8. 
D.R. Bates and M.J. Seaton show that the electron density of the D region at 
the height of 90km is about 15x10'/cc.. This is almost equal to the value derived 
here. 
§ 10. Daily Variation of Electron Density 
It is natural to consider that the electron density at night is smaller than that 
in daytime. J.E. Best et al. and A.H. Benner snow that the absorption at night is 
is about one-half or one-third of that in daytime in unit of neper. The height of 
reflection increases about 10km-15km at night than that in daytime from the observa- 
tion of two waves.. These reflection height is at the level of 92km-100km. But, from 
a point of view of the recombination and ionization, it does not appear that maximum 
electron density of the D region or D region itself exist at such a low level at night. 
H.G. Booker and L.V. Berkner®” (1938) discovered the discontinuity between the E layer 
and its lower part. They adopted, as the explanation of this discontinuity, the 
suggestion that the ion layer situates under the E layer. From this fact, it may be 
considered that the reflection of the long wave is made by the ion layer. 3 
7 
" 4 
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A study of the reflection by the ion layer will be made in the second report. 
§ ll. Variation of Height of Reflection of Wave 

In the preceding section, it was found that it is suitable to use the from 
Nye? (1—2—-e*secX) as the electron density distribution of the D region. It is interest- 
ing to study whether the variation of the height of the reflection which is easily 
derived is consistent with the measured value in the observation or not. Results of 
this study is shown in Fig. 9. Fulll line shows the measured value by the wave for 
16kc/s. and the dotted line is the calculated value. The coincidence of the two curves 


2 is very good for small secy, and it 
5 10 becomes bad as sec x increases. Never- 
= 8 theless, the difference between them is 
ie 6 much smaller than lkm. Therefore, 
> when we consider such an electron 
ES 4 distribution as to explain only the 
& wh > variation of the height of the reflection, 


Nyt (1-2 *secx) is appropriate. But 


0 0) there is a fact which we must note. 
SEC X 
Fig. 9 Variation of reflection height with secy. 
Full line represents the measured varia- a e 
tion for the wave of 16kc/s with the density distribution of Nje?t—2—© “sen, 
angle of incidence of 30°. Dotted line this curve coincide with the measured 
represents the calculated variation, for 
the wave of same frequency, using the 
distribution of _ electron density hair line. It is, therefore, rash to 
igdi—z—e-* . . - : . 4 : 
Not vs, and the hair line re- consider that the D region is the tail 
presents the calcalated variation using 
the distribution of electron density 
Nyon oo ee 2), the height of the reflection. And it is 


probable to expect that the tail of the D region has a modified Chapman’s distribution 


When we study the variation of the 
height of the reflection for the electron 


too. This is shown in Fig. 9 in the 


of the E layer, from only the result of 


of the electron density. 
§ 12. Conclusion 

It is found that the reflection of the long wave at the frequency of 16kc/s and 
170kc/s occurs at the level x=sin’O[y?,+Z-]?, and not at the level x=sin’0, or x= 
sin?@-Z. The maximum electron density of the D region is 10‘x/c.c. in daytime in 
the winter, and 2x 10'/c.c. in daytime in tde summer. The electron density distribution 
is of the type which explain the variation of the height of the reflection and 
expressed as e#(1—2 —e~*cecz) This implies that the electron density at the height of 
90km and 75km is 1.8x10!/c.c. and 6x 10/c.c. respectively in daytime in the summer 
and 9x10°/c.c. and 3x10°/c.c. in daytime in the winter. 

In conclusion, the writer wishes to express his hearty thanks of Prof. M. Hase- 
gawa and to Dr. K. Maeda for helpful suggestion during the course of this study. 
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World-Wide Character of the Progressive Change in 
the Disturbance Forces of Geomagnetic Bays* 


By Naoshi FUKUSHIMA and Hironori ONO 
(Geophysical Institute, Tokyo University) 


Abstract 


The vector-diagrams illustrating the change ‘in the disturbance 
forces during the whole course of geomagnetic bays, show the same 
characteristics everywhere in the middle latitude regions; namely the 
disturbance force vector rotates clockwise with time in the forenoon 
and counter-clockwise in the afternoon, and the above characteristic 
corresponds to broadening of: the positive bay area over the earth’s 
surface during the course of bay. It is also noticed that the average 
current system of bay disturbance tends to shift eastward on the whole 
during its development and decay processes. 

In the disturbance force of bays observed at Sitka, the influence 
of the auroral zone current is clearly seen. On the other hand, the 
disturbance force at Honolulu is so small in its magnitude that the © 
bays’ character is less systematic there compared with those‘in the 


middle latitude regions. 


§ 1. Introduction 

Vector-diagrams of the horizontal disturbing forces during the course of bay 
disturbances have been examined by L. Steiner [1], H. Hatakeyama [2] and J.Ma Princep 
Curto [3] with the aid of a number of data at O’Gyalla (47°52’/N, 18°11’E), Toyohara 
(46°58’N, 142°45/E) and Ebro (40°49’N, 0°31’/E) respectively. A well-known character of 
the mean vector-diagram, which is common for the above three stations in the middle 
latitude regions, is to form a loop during the course of the bay disturbance, and the 
sense of its rotation depends upon the local time, namely clockwise in the forenoon and 
counterclockwise in the afternoon. | 

It will be examined in this paper that whether the above-mentioned char- 
acteristic of the looping phenomenon of the bay disturbance can be also seen in 

the records of American middle latitude stations, Cheltenham (38°44’/N, 76°50/W) and 

Tucson (32°15'N, 110°50’'W). The observational results at Sitka (57°04’N, 135°20’W) 
and Honolulu’ ‘(21°19/N, 158°04/W) were also examined in order to study the character- 
, istics of the bay disturbance at the higher and lower latitude stations. 


* ‘Conitibution from Division of Geomagnetism and Geoelectricity, 
Geophysical Institute, Tokyo University. Series II, No. 34.. 
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§ 2. Data and method of analysis 

“ Magnetograms” published by U.S. Coast and Geodetic Survey were available 
in this study. The data during 1947 and 1948 were used for Cheltenham and Tucson, 
while for Sitka, the records during three years, 1946-1948, and for Honolulu, Jan. 
1946 through June 1949, were analysed. 

The whole duration of a bay disturbance was divided into six parts; the time 
of its beginning is denoted by ¢ and that of its end by é, while ¢; means the middle 
of the whole duration. The disturbance force of bays was defined approximately to 
be the deviation from the straight line connecting the values at f and z# on magne- 
tograms. Whole selected bay disturbances were divided into eight groups according 
as #,;s belong to 0-3h, 3-6h, — — -, 21-24h in local time. Numbers of selected bays are 
tabulated below. 

Number of bays during respective intervals of local time 


Observation 0-35 3-6h 6-9h 9-12h 12-154 15-185 18-21h 21-24h Total 
Cheltenham 41 42, 18 12 15 26 29 49 232 
Tucson 45 Bi. 5 11 22 31 39 48 228 
Sitka 56 50 18 12 8 10 39 80 2¢3 
Honolulu 53 20 5 38 20 40 33 69 279 


§ 3. Bay disturbances at Cheltenham and Tucson 

The mean vectors of horizontal disturbing forces of bays in each eight time 
intervals mentioned before for Cheltenham (50.1°N in geomagnetic latitude) and Tuscon 
(40.4°N in geomagnetic latitude) are given in Fig. 1. Diurnal change in the direction 


an ch h h 
0~3 3~6 6-7 9~12" pa-15 * ts2 pe 21* 24-244 
Fig. 1 Vector-diagrams of mean horizontal disturbing forces of geomagnetic bays at Cheltenham 


and east-west.) 


and Tucson. (Rectangular axes in the figure are referred to the Magnetic north-south - 


of the disturbing vector has quite the same tendency as those reported hitherto by 
several investigators; namely, it directs northward at. midnight, westward in the 


J 
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morning, southward about noon and eastward in the evening, showing counter-clock- 
wise turning in a day. The sense of the rotation of the loops in the vector-diagrams 
is clockwise in the forenoon and counter-clockwise in the afternoon at both observa- 
tories, with only an exceptional one for 0-3h at Cheltenham. Therefore, it might be 
Said that the sense of the rotation of the loops just mentioned above will be a general 
character for every station in the middle latitude regions. 
Mean angular velocity of the rotation of horizontal disturbance vectors near 
the middle of the whole duration of bays, is derived for each loop from the force 
w py vectors at f; and ¢,;, and denoted by w. 
This angular velocity w is taken positive 
Cheltenham or negative according as the rotation of 
the loop is counter-clockwise or clockwise. 
The relation between w and local time A 
(replaced by longitude) is shown in Fig. 2. 
If w assumes to be expressed by the form 


W=Q+¢,Sin(A+a,), 


the relations between w and 4 in these 
two cases are 
w=—3+17sin(A+150°) °/hour 
for Cheltenham, 
w= —3+24sin(A+169°) °/hour 
for Tucson, 
where the calculations were made under 
the assumption that the duration of the 
mean bay disturbance is two hours in 
every case. On the other hand, the 
similar expression for the data at Toyo- 
hara (36.9°N in geomagnetic latitude) 
during Aug. 1932-Aug. 1935 [21[4] shows 


w= —9+46sin(A+197°) °/hour. 


Fig. 2 Relation between w and local time > for 
EE ey In each cases, the minimum and maxi- 
mum values of w take place in the morning and in the evening respectively. Dynamic 
behaviour of the current system of bays, which is derived from above-mentioned 
results, will be discussed in § 4. 
§ 4. The bay disturbance at Sitka, and the current system of the average bay aie 
The results of similar analysis of the data at Sitka, Alaska, as well as the 
variation in Z-component are given in Fig. 3. The geomagnetic latitude of the 
observatory is 60.0°N. We notice here that the horizontal disturbing vector does not 
change its direction smoothly counter-clockwise with local time as the cases for 
Cheltenham and Tucson, but almost southward throughout the forenoon, northward 


at 12-18h and westward at 18-24h. 
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Fig. 3 Vector-diagrams ‘of mean horizontal disturbing forces of bays ie simultaneous variation 
in the vertical intensity observed at Sitka. 


Now suppose that the horizontal disturbing force of a bay at a locality on the 
earth’s surface is given by such the vector as given in Figs. 1 and 3, so far as the 
geomagnetic latitude and local time in these figures are taken as the coordinates of 
the earth’s surface to indicate the locality. We shall call a bay thus defined “ the 
average bay.” The equivalent overhead current arrow maps for three stages of the 
average bay are given in Fig. 4. The three stages are the developing stage, the 
maximum stage and the decaying one, represented by the mean state of 7, and ?é2, 
the state at ¢;, and the mean of f, and é;, respectively. Stream lines in the maps are 
SO drawn that 2/3 of observed H- and D-components. are taken as. external origin 
parts of respective components, and the simultaneous variation in Z-component is also 
taken into account. The stream lines, which divide the middle latitude regions into 
positive and negative bay areas, are drawn with thick lines, where positive or negative 
bay area means the region in which positive or negative bay is observed respectively. In 
the middle: latitude regions, threre can be seen the tendency that the positive bay 
area becomes wider and the current system shifts eastward on the whole with the 
progress of time during the course of the average, bay. This characteristic of the 
current’ system: can be explained by the relation between w and 4 in the preceding 
paragraph, and a kinematical interpretation was proposed by one of the writers [4]. 

“One can see in Fig. 4 at once that the effect of the auroral zone current is 
appeared in the disturbance force at Sitka, though. the centre of the auroral zone 
current situates higher than 60° in geomagnetic latitude. The overhead current 
system in the polar cap could not be obtained in this study. However, a current 
system such as given in Fig. 5, which illustrates the average of the current distribution 
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Fig. 4 Progressive change in the overhead cur- 
_rent arrows and the current system of 
the average bay disturbance, ~ : 
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Fig. 5 The mean equivalent current system 
of polar magnetic storms. (after T. Nagata) 

for polar magnetic disturbances [5], 
seems to well harmonize with the cur- 
rent system of the average bay in Fig. 
4. Briefly speaking, the diurnal change 
of the direction of disturbing forces of 
geomagnetic bays at Sitka seems to be 
strongly subjected to the longitudinal 
asymmetricity of the auroral zone 
current. The aspect of. the current 
system, which is obtained by combina- 
tion of Figs. 4 and 5, shows good 
agreement with the average current 
system of bays by H.C. Silsbee and 
E.H. Vestine [6]. 


§ 5. Bay disturbances at Honolulu 

As will be seen in vector-diagrams 
in Fig. 6, the magnitude of disturbance 
forces at Honolulu (21.1°N in geomag- 
netic latitude) is much small compared 
with those observed in the middle lati- 
tude regions. The disturbance force vec- 
tor directs almost northward or south 
ward, showing that bay disturbances — 
can hardly be recognized in the record 


of declination. The rotation of loops 


of vector-diagrams at Honolulu does 
not show such a systematic character 
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Fig. 6 Vector-diagrams of mean horizontal disturbing forces of geomagnetic bays at Honolulu. 


as in the middle latitude regions. This may be perhaps due to an ambiguity in defining 
the disturbance force of bays, since the disturbances themselves are small there. 


§ 6. Conclusion 

In this statistical work, it was obtained that the vector-diagrams, which illustrate 
the change in the disturbance forces during the course of bays, show the same 
characteristics everywhere in the middle latitude regions, in America, in Europe and 
also in Asia; namely the disturbance force vector rotates clockwise in the forenoon 
and counter-clockwise in the afternoon as described in § 3. The above characteristic 
corresponds to broadening of the positive bay area during the course of bay as illustrated 
in Fig. 4. It is also noticed that the average current system tend to shift eastward 
a little on the whole during its development and decay. 

The disturbance force of bays observed at Sitka is fairly affected by the 
neighbouring auroral zone current. The current system in the polar cap will be 
supposed to be a similar one to that of Fig. 5. 

The result of frequency. distribution of occurrence of bays suggests that the 
disturbance force at daytime is rather overestimated and that at night rather under- 
estimated in “the average bay.” Hence, the mean current intensity will be somewhat 
large in the dark hemisphere and weak in the sunlit one in practice, compared with 
the current distribution of “the average bay.” Nevertheless, the above conclusions on 
the progressive change in the current system for geomagnetic bays will remain without 
serious alteration. 

In concluding, the writers wish to express their hearty thanks to U.S. Coast 
and Geodetic Survey for sending us “ Magnetograms.” Their cordial thanks are also 
due to Prof. T. Nagata for his kind direction throughout this study, and to Messrs. 
R. Sato, Y. Ono and T. Takano for their assistance in the statistical work. 


(Read Apr. 27, 1952) 
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Dynamo-Theoretical Conductivity and Current in the Ionosphere. 


By Ken-ichi MAEDA 


(Electrical Communication Laboratory, 
Ministry of Telecommunications. ) 


Abstract 


Tensor expressions of anisotropic velocity and conductitivity of 
a simple ionized gas in magnetic field is applied to the ionosphere as 
a heterogeneous ionized medium. It is shown that the vertical electric 
current in the ionosphere is suppressed down to a negligible magni- 
tude, and under this condition the specific conductivities of the 
ionosphere are derived. With these conductivities the dynamo-theo- 
retical equation is constructed and solved for some special models of 
the ionosphere. The directional conductivities of the ionosphere are 
introduced and explained. The vertical motion of electrons and ions 
in the ionosphere is considered and its importance is briefly noted 
from the viewpoint of the ionospheric distribution. 


- Introduction 

The conductivity of the ionosphere is an important means to correlate physically 
the geomagnetic and ionospheric phenomena, and at the same time, so far as the 
conductivity implies the motions of charged particles in the ionosphere, it is also 
important for the study of ionospheric structure and distribution. 

Concerning the conductivity problems there are two main points to be studied. 
The one is the estimation of number densities of ions and mean collision frequencies 
(v., ¥;) of electrons and ions with neutral particles, which are not directly observable. 
The other is the analytical problem to derive necessary formulae or to solve neces- 
sary equations in order to correlate or interpret geophysical phenomena. | 

As for the former problem, the ion density was studied by Massey“ and others, 
y, and v; were theoretically formulated by Chapman and Cowling,“ and the present 
writer presented a report“) on these subjects in connection to the condudtivity problem 
and there the equation of conductivity was put simply as 


en ve €?(21 + Me) Vi 


The ionosphere, however, is not isotropic on account of the influence of the earth’s 
magnetic field. And moreover the effect of vertical drifting force recently studied by 


some investigators“ should be taken into account.. 
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The present paper concerns the latter problem. The writer formulated the 
dynamo-theoretical relation among ionospheric conductivity, earth’s magnetism and 
atmospheric. motion, taking into’ account. the effect of vertical drift and anisotropic 
conductivity of the ionosphere and then the directional conductivities, such as south- 
ward, eastward and cross conductivities, are introduced in the dynamo-theoretical 
sense. 

On the other hand the latitude distributions of thé range of diurnal variation 
of geomagnetic field have been analyzed by some geophysicists“?“, who pointed out, 
for example, an anomalous increase of the range near the magnetic equator. And 
more over it is well known that the world wide distribution of the ionosphere conforms 
with geomagnetic coordinates rather than geographical ones. 

The above-mentioned dynamo-theoretical formulae with directional conductivities 
are solved for some model cases as applied to E and F layers and the -ionospheric 
current density calculated from the above solution shows a possibility to interpret the 
geomagnetic result of analysis above-mentioned more rigorously than by the usual 
dynamo-theory. 

The vertical motions of electrons and ions derived from our present theory are 
found to be significant in their effects on the density and world wide distribution of 
the ionosphere. 


If. Velocities and Conductivities of an ionized gas in magnetic field 


As the problem stated in the above title was treated already by some investi- 
gators®®, the writer will review briefly the problem. 
The equation of motion of a charged particle with a mass m and a charge e, 


during a period between its successive collisions with neutral particles, is as follows. 
Els: | oO = E+[v-0n] 
i 2" m mir 


where v is the velocity, E is the external electric field and ,, is a vector expression 
of gyro-magnetic angular frequency. When His the total magnetic field, w,, is equal 
to eH/m. 

When no external forces act on the charged particle, the total sum of the instanta- 
neous velocities of all charged particles, Sv, can be assumed as zero, or the average veloci- 
ty, v( = ne v }, is zero, m being the number density of charged particles and ¥ operating 
to a unit Pounce When some external forces such as the electric and gyro- magnetic 
fields are present, the persistence of velocity after collision should be considered. 
If the free electron is the charged particle under consideration, the persistence can 
be neglected and v of an electron at the instant immediately after collision as 
averaged over all possible magnitudes and: ‘directions of the neutral particles collided 
by the electron is zero and therefore we can put Sv(=nv) as zero. In the case of 
negative or positive ions, however, the persistence ratio (w) of velocity is not zero, 
but ranges from 1/3 to 1/2 according to Chapman and Cowling®). In a rough 
calculation it can be allowed that Sv(=nv). is solved first by putting Sv(=xy) after 


<< 
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collision as zero and then the part of velocity increase accumulated at each collision 


is added to Sv, that is, approximately (1+2)Sv. Q is the sum of a series of w and 
roughly from 1/2 to 1 (near 1/2). (1/2 and 1 correspond to zero~ and infinite 
velocities respectively.) Therefore the equation (1) can be simply thought to give the 
velocity of forced motion of a charged particle with zero initial velocity as average. 
Then the solution of (1) will be written as v(t, to), which is the velocity at time 7, 
after the charged particle undergoes the last collision at time é, If we put the mean 
collision frequency as v, Sv is given as follows. 


af = t 
C2) SU=nv= n| U(E, f))e-*¢-*o vat. 
“ppt v(t, to) is solved from (1) for the follow- 
2 


ing general case as shown in Fig. 1. The 
notations are defined in the followings. 


x,y, 2 : three axes of rectangular coordi- 


magnetic 2 

field nates, and in the later sections 
eastward ot 
% (18nd they are taken respectively as 
southward, eastward and upward 
x - mass velocity directions, that is, 70, rsin@-4 and 
sodomrd z, where 7, @ and 4 are radius, 
Fig. 1 Coordinates, mass velocity and colatitude and longitude of the 
magnetic field. point. 


X, Y, Z : three components of magnetic field, 

U, V, W : three components of mas velocity, (V), 

E,, E,, E, : three components of electric field. 
The last six quantities listed above are assumed to vary slowly or to be nearly 
constant with regard to time. 

0, Vy V, : three components of velocity of forced motion of a charged particle 

by electric: and gyro-magnetic forces (peculiar velocity). 

@,,Wy,0, : three components of wo, 
(e, denoting the electric charge, contains the sign, that is, w’s are negative for elecizons 
and negative ions and positive for positive ions.) 

With these notations the equation (1) is solved and three components (Vz, Dy» V v.) 
of average velocity v are obtained as follows. 


—- @ vito,” “ | vo,t+@y S , —vo,+o,0, A 

Se RO Se PNA oR PORE Fale 

—  @ f{—va.+ 00,4, Y+o,7 4 YO, + OO, B } (i 
ot) a { y+toz a y? Tpak? yt ton, "J’ 

— YOy+o,0, A gp ct Oye V+o, rat 

if mv { +o; <4 ey See yr Pon — 

ax B =E,+ZV—YW, : 

63’) Sealy By=Ey+XW—ZU, 


E. =E,+ YU-XV. 


current density by (7), where the summation symbol + means the sum of the thre 
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Then the conductivity is obtained as follows. When é-z is multiplied to both 
sides of (3), the left hand side becomes the current ty forced motion and its three 


components 7,, 7, and 7, are obtained. 


A 
i,=CNU,= Oxe* 7s + Ory" Ey + Ox2 “E. 
OAD ee bis Oe Ove" Ee + Oyy* B, + Cys" E., 
|, =eNv,= One" EB + Oxy" B+ O22" BE. 

where Ox2 Ory. etc. are the elements of tensor of the specific conductivity of an 


anisotropic ionized gas. 


en +0," _ @n Vo,+0,Wy re Cn —VO,t+O,0r 
Oxx = a5 2° 40 2 2 ’ a poe 2 ’ 
MY V+Wm my v+Om my y-+Om 
(5) - en —VvoO,+0,0, eer en v*+u,? A CN VO,t+W,0z 
"=" m Veo * "my v+o2’ ° my Pe 
oe CN VO y+O.0x _ On —VWz+ Wy, ae en v+o,? 
= my Po a pi gk? OF my eb oee 
The velocity and conductivity for ions are to be taken as about 50% greater than those 
given by (3) and (5) respectively according to the consideration of velocity persistence 4 


mentioned above, but as the numerical values of x, »y,... etc. are fairly rough, this 


factor (5024) is neglected in the later calculations. 


III. Vertical current in the ionosphere 

The ionosphere contains besides neutral particles the three kinds of changed 
particles, that is, electrons, negative and positive ions, whose number densities are 7, 
m, and m and masses are m and M (common for both ions) respectively. Then “- 


the volume density o of electric charge is given by (6) and the three components of Te 


__ kinds of charged particles. . 


o=V @ (m—,—2); 
=2%,+0U, L= PU ease ole =21, +oW. 


=q-ann— —Bn +rm,—div (ny), 


4 xia 


i} hs d 
Aan oy Sn a (am), sp m 
= i« ees ; 


Se 
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obtain the following equation of continuity. 


‘eg 0 healt ol LTA 
ie ak Oran Ady Gens 


And moreover Poisson’s equation must hold as an additional condition in the potential 
field with the presence of p. 


1 (9E, , OE, , OE 
(10) a ( 24 OEy 2 
ae i ae ee 


Then from (9) and (10) we can derive 


anb (Gt) - oo: = dnl + al (F+ : Gyr) a2 +{{ S. ( ie - (2) laa 


and then replacing J, of (11) by the third equations of (7) and (4), 


(12) Ee Sasser: {erie 4 Sean P: 
dat ins) 225 
where 
A A 
(13) I= D316 22° E+ Oy: Eyt+d3022°(YU-XV)+Y 2(n2—m—n) W, 
0 (dE 0 (dE,\)\ OL af, 
14 ==) : y 
oe a Hiselcee aaa Nae +(e dy ae 


dE, 


When we assume }}0.,, as constant and J, P and ( as slowly varying with 
- : . 2=0" . 
a time factor e7**, where w is 0.728x10-4(sec~1) (24 hour period) or two or three times 


of it, we get as a solution of (12), 


4n(§-+ P)— toe iat 2 | 
Cacti TS eee: git __g-4nBo , 
Gay ae aeichd 


4n>1022__f nd Tees 
yp ed ALE eee p= 4nEo st 
Msg do | ( dred cule ss 


JOJG—An>d}6,2-P— -io(@) 
at z=0g jet 
An )\022+jw 


z 


In the case of ionosphere the time constant iz eer is very small (say 10-° sec., but 
this value is not valid for the estimation of time of decay of the transient term, 
because d,, etc. were derived on the assumption that the time variation of phenomena 
is very slow in comparison to y or w,,.) Anyway as we discuss the steady state 
phenomena, the transient terms can be neglected. Then we have 


(16) cae: Saat 


where w is neglected in comparison to 47>}d,,. 
In order to estimate JI,, P given by (14) is considered. The part of displacement 
current dE, dky in (14) is very small in comparison to the last term of (14) 


dt’ dt 
(roughly 10-2 times when w=0.728 x 10-‘sec~'), and therefore P can be put approxima- 


tely as follows. 
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OF set Oly 
= od dz. 
ae | ( Ox o) 
This equation can be obtained also from (11) by simply neglecting the terms of 


displacement currents. Thus we have as the order of L, 


OU) ~O(P)~0( 1.) ~O(441,) 


As r is the radius of the point under consideration and 2, is the equivalent thickness 
of current sheet, the order of J, is about 10-!*emu, when J (current function)~6 x 
10%emu or J,, Z,~10-'emu. Therefore the vertical current is checked down to about 
one hundredth of the horizontal current, that is, J, can be assumed zero in’ the first 
order approximation and hence 2. is practically zero as well. Then the equations 


(16) become approximately as follows. 


A 
we 55 be Sees —~(YU—XV), : 
ae’) or fe Eee Bet Da By 
z : 5, D530 2s ‘ ¥ F q 
pe =a, 1 


where P is cn in Boinoaicod, to J and ye gestae et is neglected as exptane 
ed in the next section. (Ge — +, 


i ~— 5 
i ; 


Vizyi Dynamo-theoretical relations working in the ionosphere 


As mentioned above the time variation of electric charge density is zero. To © 


estimate the order of 0 we employ eq. (10). The order of £,, Ey is roughly estimated _ 
as 10-° esu, using the values in the foregoing secti d or 

ee going section and 10? esu as the order ee a 
and that of aa By ae 104 esu. The order of E, is about 10-” esu from (15), 


{ 
8. 
éx and that of’ a can be conan as less than about 10- 18 (maximum of | ; 
k= put as 10-” esu/10km= 10" esu). Therefore’ ‘the order of o is about 10-8 esu . tha 
=" aad zero as bi a with the charge seme of each Sieh a abies (en, 
<e : ~ 
= oo 


ye L=0o eet Tart Ew 
I, sca s O' yy 
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vertical current is taken into account. 0’,,, 0! zy-.etc. may be called the elements 
of tensor of the modified anisotropic conductivity of the ionosphere. 

| In the usual dynama- theory the.current function (J) is introduced .to connect 
the Sq field with the velocity potential (v). The conception of J contains the assum- 
ption to treat the ionosphere as a: single thin shell. But the ionosphere is to be 
conceived as an accumulation of many elementary thin shells, the properties of which 
depend upon the height. Therefore in order to discuss the detailed behaviours of 
charged particles and to catch integrated properties of their effects, it may be 
convenient to introduce a unit shell (or unit depth) current function j j Meaaagaete ees ts 
g [, and I,. j is defined by the ie ; 


mas ens lor J=\ jae, 
i TE Oz : 
; 1 69 -— tet 405, 
con aind0, fot td 00'S | 


As I, is everywhere zero in the first order approximation, there exist vite flow in and 
out of electric current in the vertical, direction. j is solved for each elementary thin 
shell and can be employed to discuss behaviours at the specific height of the ionosphere 
“St a J can be. obtained by integra dna. over the entire height range of the ionos- 
phere. 


i 


| » £, and E, are connected with a static potential S as follows. 


Ee= Gh. me: : 


j i 3 a “En 4 as ry A OS. | 
| : 


Z Pe 
( pee Se ae: \" i 
} } f i. 


When we sake Fe axes x, y. bore z as southward 0), ‘eastward. ine d upward 


(7), we’ obtain ‘from he the following fo nulae by liminating SE i Roe 

TE cs) Bs )j ‘J. Oo (ee OF rT bl Oo ae bbe OF si 2b ay Spal 
fe oe Ta ) “al *+ OO owe 0 ae OG OS | 

cnt. ae aes KT a fers eee 


Bei ws Fy : Moree chd os dy Pe Oiler pd 
= or ts ph Zi ole Savvy d taint Wind), 


Beer ot an Perricadk foe 


70 K. Magpa 


(21) 
O,Wr YOr 0, \2 

v+o, { ee | vo; UG: 2 z UO a = 
Ons=DI04 P+oe “oS +o? + Dd, 2 aes 0 +o? 
U6 rao? +o 

vo, )? 

Safe ts ons [27srod| 

"y+ 08 Data ; 


Own = D)}04 z 2 2 2 xt 2 
V+ Wm Slo.— Tos Gis No rad 
vto, P+ wm 
OO) 
a} 22 
yo v+W, ee {2 ers 
"y24 w Mo v+tw,? 
"y+ we 
Dong, NO, * 
Vw, v+o2 
O09 = — 00 = 210, = rer a vio? 
| Bot ok 


O,0, 
vito? | Danee aay 
+) Bora or +o? 


Ort et 


o.0 
oO zz 
Ong =} 53: VETO ahaa 


en 
the Jc itudinal ad and = ele fo ae 
where d, is longi CO} uctivity or electrons, 27 “Me for negeted 


pee ne oa Hot positive ions and v's mae w’s should also ‘be taken each value 


responding to the kind of charged particles. ‘Then the dymamo-theorstical ‘mul: 
48) becomes as follows . CARS ae 
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identical to respectively 6',., O'yy... etc. in (17) (modified conductivity tensor). As 
is shown later, the terms of d,s and de, in (19) and (22) vanish, and dyg and dyr 
for F-layer are practically equal to o’,, and O'y, respectively. For E-layer 0'zy 
and o’,, are not generally zero. The terms of dy 9 and de, however, can be neg- 
lected as an approximations, as indicated in the next section, then Oys and Oywx 
correspond roughly to o’,, and 6’,, respectively, but they are not equal to each 
other. Here we call dys and dy, as southward and eastward conductivities respec- 
tively. 

The permanent magnetism of the earth is assumed as given by a dipole magnet 


with its axis coincident with the earth’s rotation. Then follows the equations, 
a=Gesin’?, Y=0, Z=2Gcos0, 
(23) H (total field) =GW1+3cos*0, 
G= es emu. 
3 
Some numerical calcula- 


tions of o's for E and F ! 


layers were carried out. 


The numerical values of @ 


necessary quantities em- 


@ Sen-Ge) a4, Y=1210% 
i. =r 


ployed in calculation are 
given in the following table. 


ll 


at 
Ma 


Ge) K=14 % =/x10?, 
=a 


20 10 4) 


90 80 70 60 50 40 30 
colatitude in degrees 
Fig. 2-A o’s of E-layer. 


Table. 1 Elements of model ionosphere. 


on 28: x 
layer | ” ee Ye Vg yan (14+2y) m/M | ome | ! | ome | 
- = SSS eee eens a eS BE 
‘ 1x105 | 2108 6.5 | 14 
® i 0.2 LA 
E | 10° 34.5 | 70 
OA canes ees pales ae Bia" ra X10! 6x 1081 +3 cos? | 2 x 10%)1+3cos? 6 
ee ar = pee tp 20 
fp A @ 1x103 25 0.5 2 
oe | 20 
Qg) 2x10? 5 (ix 2 
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The calculated results are 
shown in Figs. 2 and3. The 


comparison of Oywx and d'y, 


for E layer is also given in 
Fig. 4. 


V. _ Definition of model 


== 
a z 
Ls 


23 0. 


ionosphere 


As seen from _ the 


en---- above figures in the forego- 


ing section the following 


points are remarked. 


(1) For F-layer dys is much 


greater than Oy, near 


the equator and there- 


fore the term of Oys 


in (22) is negligible as 


& 


ae 


compared with that of 


Ls) 


Gs in esu. 


Owz near the equator. 
(2) For F-layer 0’2, and 0'y, 
in (17) are practically 


wa Se 


ow a Sa iS) 


i 
Pig 


zero, so far aS we can 


assume v,;2<w,,; and »v,? 


<wm, and | O20 | (| Gaal) 
in (22) is very great. As 


the consequence the 


terms of dy. and de, 
vanish for the entire 
range of @ and A. 

(3) For E-layer dys is fairly 


16 greater than O,y,z for 
colatitude in degrees the zone of low latitude 
Fig. 2-B o’s ‘of F-layer. and therefore the term 


of dys in (22) is negligible as compared with that of dy, for low latitude zone. 
(4) For E-layer | dx | (| Ge, |) is comparable with o,, except near the equator and 
except for greater values of «, but except near the equator it does not vary too 
much for @ and 4, and we treat it as independent of @ and 4 as an approximation, 
then the terms of dy9 and do in (22) vanish. Near the equator | dx» | is very 
great and its term in (22) can be neglected. Therefore for the sake of simplicity 
we neglect the terms of 69 and dg, in (22) for the entire range of 6 and 2. | 
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N=10°, KeLg, 
Memo? 

4, = 2110 /+3cos'O, 
= 6:0 + 3c0s, 


——— 


Sy and Gye in esU 


Ccolatitude in degrees 


Fig. 4 Comparison of 0’,, and Oy,» of 
E-layer. 


Hee 
eo = ga 70 60 50 40 30 20 10 0 
colatifude in degrees 


Fig. 3 Oys/Owyz of E and F-layers. 


In accordance with the above- 
mentioned points we define our ionos- 


phere model as shown in Fig. 5 and 180° 
below. 
ame ieoas trom 0=0° to 0,;= 15"; 
Zone Il -as from 0,=75° to 0)’ =105°, 
Zone Ill as from 6,/=105° to @=180°. 
For zones I and IU, 

F-layer : > Oys=Owzr=Oor, 

E-layer : Oys=Owz=Oor; 


Ocr and dog are average values ap- 2x (Gus = Swe = Foe ot Soe 
propiate for the zones of each layer. @& 
rs) E-layer 

For zone II, s 

F-layer : Owz=500r Ons> Owe: » 

E-layer : Oye s00m Ovs>Owe- 
As dys OF Owx is proportional to the 
density of electron, the consideration » Ogg~ga a 50 330-2010 0 
of latitude distribution of monthly  colatitude in degrees 
average of n{(f°F2) or (f°E)?} is Fig. 5 Zones and o’s of E and F-layers. 


necessaty. According to H. Ueyda“™ eee curves are rouph oe 


tions of Oyx of E and F-layers. 
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the geomagnetic latitude distribution of fF, averaged over one year (1945-1946) is 
shown as a ratio (k) to the equatorial value and the calculated #° is also shown. in 
Fig. 6-A. In Fig. 6-B dwz, when the electron density of F-layer (#) is taken as 


40 ih 

0.9 ss) 

0.7 7 

> 0.6 (esu 6 
# a5 Ls 
0.4 2 ° ° 2 4 

Bs f= {§;/#;,(0-20} ; 


res 80 70 60 50 40 30 20 40 O 


Co/atifude (degrees) 

Fig. 6-A Latitude distribution of one year Fig. 6-B 
average of daily mean of /f°Fe 
(1945-1946). - @) 

k: ratio of f°F2 to that of equator, (2) 
Rk: calculated from k. 


independent of 6 (10°cm-"), and dyz, when x is put as 10° x ecm, are shown, saci 3 


for the case where | k= we y,=103 Bets ete 


follows, putting vertical mass motion (W) as zero. 


For zones I and III of E and F layers, 


ee WRN SEE $x atin OF. eae 1 4 
: Tr ’ Gox\Son) OA sin 004 7 Fox For) 


at pS aes th 


°95-B0 70 


S-(n=10 em) 


60 50 40 30 20 10 
Colatrlude (degrees) . 
Latitude distribution of ¢yz of — 
F-layer. tee 
n is assumed as independent of 6, 
mn is taken as pase oe of 8, ity 
108-#2cm-3,_ 

i is given in the: left figure.) 


Pak e. 


yz shay 6) fr Sz aah 


ee a 
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hand side of (28) or =. is calculated. 


(31) Anish (ZV) +—7 oy (ZUsind)\ = 2Gsin 0S} a Pm +1(9) cos (mt + om). 


m= 


Next we put the unit shell current function j as follows, 
for zones I and III (E and F layers) 


——? 
(32) ai Sox) 33K pia 


Jir=—2G(0or or Fon) IRR 


Pm+1(0) +4,t,(0)|eos (mt+omn), 


Pr 1 (0) at. a,0n(0) cos (mt ae ads 


(m+ eee 
Se'sin® 0 

(33) Hm O= TF cosey” 
22S, 


Ym (== cos ay 
for zone II (E and F layers) 


(34) Ju=—2G(3d0ox" or 500 ame L 


0 0 ‘ms 
Sm 1 OnlO+ Duh (O) + dq cos mt +08) 


._The functions 2,,(), »,(0) and h(0) are introduced in order to satisfy the boundary 


conditions between zones II andI or II. ,,(0@) and v,,(0) were first given by Rikitake@” 
in a form somewnat different from (33) as a solution of (28) with its right hand side 


zero. h(@) is a solution of (29) with its right hand side zero. Qn(9) is a particular 


solution of the following equation reduced from (29). om ( 


(a. OMe bee ea 
a0 $y {sina Se pes (0). 
@,(0) and h(@) are given below. — ia er 
| 2 m421 (mes sd 
H ——— eee u 0, ‘ 
ey Qn(0)= Wieser gmt. ce ae Cet 
! co “a 
(36) : _ — hO)= lon tes. ee. 


The constants Qn. b,, and d,, are calculated from the boundary conditions listed below. 
‘ a ee 
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Table 2-B. 
m(m+2) 
-n eee Om(9) | dy 
Oe 
sinOcosi 8 
t 2 ae: | n/4 
2 ee sin”6 cos 8 + — cos 6 | 0 
ey . 9,410 sindcos®  9)/10 98. 9/10 = 
3 sin acos 6 + 16 pit € Luar 16 5) 16 Ti 
pad te By 55. Je 41/35 
4 een 8cos 8 + 5 sin?6 cos 9 + 15 cos 8 = 
Numerical calculations were carried out for m=1 and a, and dy, are as follows. = 
For F-layer : ~ 
ay—0.152; 5b, =0.716, 
For E-layer: Fk 
a,=0.10, 5,=0.635. | 4 
¥y 


The current densities (I, and I,) of southward and eastward directions are calculated 
with the relations to the current density function (j).— 


Therefore for F-layer, 


"1 7 = 2CSorkas! cost + 0.1522 |sin (¢-+.432), . sini’ <4 
r 1l+cos@ J - < 

(37) aay Br te . forvzone 5: bh) suliee 

7 ork {COS 0.152 Y : ie 
¥ T ae te 1+cos@ alee 5 
Tt 
is 10G donk! {cos 6 05(5-0 ) 0.716, 
a) + (a oslo! eee —— logtg— sin (t+,5), 

(38) PHG9 2 sin 0 sind 2 

ger - : 
a 10GdO ork; es 20 manne 0.716 cos (t aye’ 

Yr 2 sind Jf - 


and for E-layer. stares * 
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Bia +0. LTpapea fined, 


1+cos@ 
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— 6,500 km, Oor= 10" esu. 


Fig. 7 shows the latitude distribution of current desity for E and F layers and there 


Oor of E-layer was taken as 10 times of Jor Of F-layer. This assumption is derived 


from the comparison of Own of each 


=-Ys ems, layer, based on the observed data of 


-2.9x10" emu, 


electron density. 

The current density as _ the 
function of height (z) is thus given by 
eqs. (37), (38), (39) and (40), where 
Oor, OY Ooy varies with height. (7 


does not vary too much.) The curves 


shown in Fig. 7 are the values of J, 
and J, at the height of their maximum. 
The integrals, \Zedz and {L,dz, from 
the lower boundary of F-layer up to 


colatitude in degrees 
° 


Iy(E)+ ze Ag) 


Z 
Za(F) _ 
Za(E) ¥3) 


the upper boundary of F-layer give 


the total current (current/length), and 
\Z.dz and (J,dz divided by the equiva- 


ry 


lent thickness, (z,), (Zedz/2q and | I,dz 


/2a, are the current density of the 


Aiurnal range of. 
70 |-Aorizontal tela 
by M. Hasegawa 
80 Harbitrary scale) 


210 08-06 -04 -02 0 02 OK 06 08 10 /2 Lh 16 x(3.3x10%) 

2 ee overhead current sheet, where the 

current densify in emu. 

Fig. 7 Current densities of E and F-layers. ionosphere as a whole is assumed as 
ig re, Ge: E-tayers 5 a single shell with z, as its thickness. 
Za(F)/Z2(E) =5/3) 

As a rough treatment the sum 

of the current densities of both layers is also shwn in Fig. 7 and this quantity is 

roughly proportional to the first harmonic of the range of diurnal variation of geo- 
magnetic horizontal component. For comparison the curves of the last quantity 
analyzed by prof. Hasegawa‘” are shown in the same figure. We find thaf both curves 

(theoretical and observed) are in good agreement. 

VII. Notes on the importance of vertical motion of electrons 


In sections III and IV it was clarified that p and ae can be assumed as prac- 
tically zero. But as for each kind of charged particles, a vertical flow does exist and 
the velocity is given by (3). To calculate the peculiar velocity given by (3) the electric 
and gyro-magnetic forces are derived from (16’), (17) and (37), (38), (39) and (40). 


E,+ZV=—2_I,9,(0)sin(¢-+ a) -—2— Iy9,(9) cos (¢+), 


Ons OX) ' 
: 0, 
E,—-ZU=— aa I,0,(0) sin (t+a) + Fz boPv(0) 008 t+ 0), 


(41 : 1 pe oseweiat 
4 ts yu-xV={ a0 Be |.) sin(t+a) 


Ons vege Oox Oz 
ke Th MeL 22 Mi 0) cos (+a). 
+{ Ose £DIGt. | Owe Don f WG ( 
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In the calculation of the above equations, J, and f, given by (37), (38), (39) and (40) 


were written as follows. 


( I,=dodo0,(0)sin (+a), 
(42) I, = Gold (0) cos (t+ 4), 
I, =2Gh,}/7, Oo=Oo0r OY 300 and Oop OF 50op 
These results of (41) are introduced into (3) and each component of the velocity is 
calculated and given as follows, where the suffix 7 means the kind of charged particles ; 


n, is n for electron, 2, for negative ion, and 2, for positive ion. 
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Especially for F-layer (43) can be made somewhat simpler, as an approximate a 
sion. 


Cire | aah. 


res do L (0) Onxt— O 224 Von 0 Id 6 Y . 
Vxi = se Ser cn ae “om of dcos(t-+-a)-Z¥* & shh 
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expressed by div (zy), diy (v1) and div (#oy=)in (8), that is as follows. 
Dg 0) _ 0 _ 
Og Mz Gz (mrz1), Gz (z2)- 


The third equation of (44) is considered for estimating v,, of F layer. For 
example — By (uve), that is the case of electron, is formulated as follows except near 
the poles. 

se 


0 _ 1 rf) Osx—Ozz a eee il 0 
45 ——=;= 70.) = = =s Didier 0: 
(4) Oz (nv:) VY e {Z. a ae La a, oh 0 “ \, 


0 : 
=AS Bee MP Osin tas!) +B Sf En \ go) cos(t+-0, 


where A and B are calculable constants, (9) and g(0) are functions of 0, and ky} /r 
has a dimension of velocity (negative value) conneted to the velocity potential. v,; is 
taken as 25 at the height of maximum electron density. 

It can be proved that A is about one tenth of B, which is of the order of 
unity, and therefore we take into consideration the second term of the right hand 
side of (45) as an approximate treatment. Then using (37) and (38), we have for zone 
i; 


Oe BO ReiN\:. cos 20 2 
46 meek, yoo "CT , & Yano Dung. eel “i 
(46) ial (x )sin (See area srereerey. eo Gas 


and for zone II, 


01S 15 0 Ry cos 20 0.716 
7 ———(nv, fobs 0 —0.5 
ie Pets 4 02 ("5 )sin 2 ae sin 0 


Dz teos (¢+4a,1). 


For ions, ” should be replaced by 2, or m, in the right hand sides of (46) and (47), 
that means v,=v,1=0V-2. 

The vertical motion of electrons and its effect on ionosphere were treated 
already by Martyn™ and other investigators“, but the writer will describe briefly 
some extensions of his present theory. 

The value of oy nes) amounts to a considerable magnitude, which can not be 
neglected in comparison to q, a”, Bn in (8). The first equation of (8) can be written 
as follows, when at a certain time of day (say noon) and a certain colatitude (say near 
the equator), 40.) is put as S(ne), and 7 (detachment) is neglected, 


3) er (an +B) np (ne). 


To make the problem simpler, it is assumed that ¢ is constant, q is Chapman’s model 
of ion production and om,+8 is put simply as 7, which is given below. 
2Zm—2 


Ota + B= = Hoe eyes 


where z,, is the height of q,,(0), that is the maximum ion pronactios at %=0) Then 
(48) can be written as follows. rips 
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a a 
(49) or +c ae ~q—yn. 


The partial differential equation (49) can be ey ee and numerically calculated, 
but here we consider the steady state of where aes is neglected. Then the differ- 


ential equation, 
=) ‘a ee 


is numerically treated with c= +1000cm/sec and +100cm/sec, H=30km, 7% =10~‘sec™* 
and ~=0. The results 

tere Ny of calculation are shown 

P as in Fig. 8. When c=0, 


et 


M is equal to q,,(0)/7o- 
The curves are drawn 


mm 


22-2 
for n/m versus we 


As seen from the figure, 
both upward (v,>0, c>0) 
and downward (v,<0, 7 
c<0) motions of elec- 
trons make the maxi- 


mum of 2 decrease. The 


maximum value of | ¢ | 


a 
4 


is comparable to | k,1/r|, 


RYn. 
Fig. 8 Vertical distributions of electron density (F-layer). that15) nearty 1 
c= 10° & 10%: upward drift of electrons, or less, with | R,' | =3x 
(- —10° & —10: downward drift of electrons, 10" emu. Therefore, for _ 
c=0 : no drift of electrons. as 


‘ | example, the maximum 
electron density around noon at the equatorial zone is about half of 2) and the height — 
is raised as ¢ or v, has a positive sign, that is, electrons move upwards around noon 
A for the above zone. ; ass piece 
Fs _ Although the above treatments and descriptions are rough, it may be 
t the diurnal variation of electron density and height of the ionosphere 


> 
nat 


distribution of the ionospheric density and height will be distorted, so as it , 
magnetic coordinates rather than geographical ones. These probl 
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return to the first deduction of specific conductivity in Section II and to treat the 
problem without the condition of slow phenomena. The problem is not simple, but 
it is supposed that the mean collision frequency (v) is the controlling factor and the 
time constant of the transient of vertical current decay described in Section TI will 
be greater than 10-8 sec, but is much less than 1sec. 

(2) The vertical current in the steady state is practically zero, as the ionos 
pheric thickness is very small in comparison to the radial distance from the earth’s 
centre. The ionospheric current is generally composed of electrons and ions (negative 
and positive), and therefore the zero vertical current does not imply necessarily no 
vertical motions of each individual charged particle. Electrons and ions’ move 
vertically with their own velocities, each of which is given with the electro-motive 
force connected with the dynamo current under the condition of zero vertical current. 
These vertical motions of charged particles play an important role on the electron 
density and height of the ionosphere. 

(3) The unit depth current function (¥) is introduced by the writer in Section 
IV. As defined in Section IV and described in the last part of Section VI, the density 
of dyamo current at any height of the ionosphere is given by 7 and the usual current 
function (J) is exactly derived by integrating ;. 

Also in order to discuss the effects of motions of charged parcels on the 
electron desity of the ionosphere, the necessary velocity of each charged particle at 
every height is obtained through the function 7. 

In this report we described the cases where the conditions v.2<wZ., v2 >w2,; 
(for main E-layer) and v,2<w2,, v°<w;, (for main F-layer) hold. We did not calculate 
the case, where v; is comparable to w,,;, that is to be applied from the upper part of 
E-layer to the lower part of F-layer. The investigation of the above intermediate 
region will be necessary, in order to obtain the current function (J), through which 


‘the ionosphere and geomagnetism are physically connected, and also in order to 


investigate the vertical distribution of the ionospheric density, its time variation and 
world wide distribution. © 
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